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ABSTRACT 

Water-rich super-Earth exoplanets are expected to be common. We explore the effect of late giant 
impacts on the final bulk abundance of water in such planets. We present the results from smoothed 
particle hydrodynamics simulations of impacts between differentiated water(ice)-rock planets with 
masses between 0.5 and 5 M® and projectile to target mass ratios from 1:1 to 1:4. We find that giant 
impacts between bodies of similar composition never decrease the bulk density of the target planet. 
If the commonly assumed maximum water fraction of 75wt% for bodies forming beyond the snow 
line is correct, giant impacts between similar composition bodies cannot serve as a mechanism for 
increasing the water fraction. Target planets either accrete materials in the same proportion, leaving 
the water fraction unchanged, or lose material from the water mantle, decreasing the water fraction. 
The criteria for catastrophic disruption of water-rock planets are similar to those found in previous 
work on super-Earths of terrestrial composition. Changes in bulk composition for giant impacts onto 
differentiated bodies of any composition (water-rock or rock-iron) are described by the same equations. 
These general laws can be incorporated into future N-body calculations of planet formation to track 
changes in composition from giant impacts. 
Subject headings: planetary systems — planets and satellites: formation 



1. INTRODUCTION 

Super-Earths, massive terrestrial exoplanets in the 
regime < 10 M®, are expected to come in a diversity 
of bulk compositions. Both planet formation theory and 
comparative planetology stand to benefit from distin- 
guishing super-Earths of different composition. It is es- 
pecially valuable to know how water-rich such planets 
could be. Water is abundant in protoplanetary disks, 
but it is also very volatile, so the details of the planet 
formation (and post-formation) process can be crucial 
in determining the distribution of water-rich planets by 
bulk composition and orbit. 

The catalog of observed extrasolar planets now in- 
cludes more than 400 membersQ Of these, more than 
70 planets have been observed transiting their parent 
star. Until recently, the only known transiting plan- 
ets were similar i n mass to either Jupiter or Neptune . 
With CoRoT- 7b (ILeeer et alJ l200l IQueloz et all 120091) 
and GJ1214b (|Charbonneau et al.l I2009T ). we now have 
the first observed transits of planets in the super-Earth 
mass regime. Many more detections ar e expected in 
the n ear future from the Kepler mission (Borucki et al. 
2009). A transit observation provides a determination 
of the planet's radius as well as the inclination angle of 
the planet's orbit with respect to the line of sight. When 
combined with the mass determined from radial velocity 
measurements, the mean density of the planet is derived. 

Given the bulk density of a super-Earth, it is possi- 
ble to estimate a range of possible internal compositions, 
from rocky/iron terrestrial planets with small radii for a 
given mas s, to icy and gas rich planets with much larger 
radii fe.g..lValencia et al.|[200l[2007l:lFortnev et all2007t 
iSeager et al.l 120071) . The precise composition of the 
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planet cannot be determined from internal structure 
models alone because of the degeneracy in the mass- 
radius relationship for various compositions. However, 
with precise measurements of radius and mass (uncer- 
tainties < 5% and 10%, respectively), one can distinguish 
an icy/rocky su per-Earth from a rocky/iron super-Earth 
([Valencia et al.ll2007t) . 

The mass and radius measurements for CoRoT-7b, 4.8 
± 0.8M® and 1.68 ± 0.09i? ffi , respectively, are consis- 
tent with a planet of terr estrial composition (|Leger et al.l 
I2009t IQueloz et al.ll2009h . GJ1214b, with a similar mass 
of 6.55 ± 0.98M e , has a radius of 2.678 ± 0.13J?®, im- 
plying the presence of a thin hydrogen/helium envelope 
around a core composed primarily of ice or, alternatively, 
a thick hydrogen/helium envelope around a rocky core 
(jCharbonneau et al.ll2009l ). 

Further constraints on the compositions of transiting 
super-Earths can be obtained fr om a theore t ical u nder- 
standing of pla n et fo rmation. Ilda fc Lml (|2004f) and 
IMordasini et al.l (|2009|) have applied detailed models of 
planet formation to t h e gen eration of synthetic popula- 
tions. iFigueira et alJ (120091) have c o mbin ed the planet 
formation model of Mordas ini et al.l (|2009T ) with simpli- 
fied interior structure models to determine a statistical 
distribution of possible compositions for the transiting 
hot Neptune GJ436b, excluding the possibilities that 
the planet is either an envelope-free water planet or a 
rocky planet with a hydrogen/helium envelope. Such 
constraints represent a major triumph for planet forma- 
tion models. However, all current models neglect the po- 
tentially important effects of multiple embryo formation 
and giant impacts. 

Giant impacts are thought to have played a major role 
in shaping the planets in our solar system. Such events 
have been shown to explain the high bulk density of Mer- 



cury in the inner solar system ()Benz et al.lll988l 120071) . 
the formation of th e Haumea collisional family in the 
outer solar system (|Brown et al.l 120071 : iLeinhardt et al.l 
I2010T ). and properties of several planets in between. In 
fact, giant impacts are an inevitable consequence of the 
core ac cretion model of plan et formation (e.g.. lWetherilll 
I1994I) . iMarcus et al.l (|2009T ) developed scaling relations 
to determine the outcome of giant impacts between 
planets of terrestrial composition They extended the 
catast rophic disruption criteria of iStewart fc Leinhardtl 
(2009) into the super-Earth mass regime and provided 
a power-law relationship for the iron mass fraction of 
the largest impact remnant given target mass, impactor 
mass, and impact velocity. These laws were then com- 
bined with simple dynamical arguments to constrain the 
minimum radius of super-Earths as a function of mass 
([Marcus et al.ll2010Q . 

Ultimately, we want to know what the distribution of 
super-Earths of different bulk compositio n looks like in 
the m ass-radius diagram. The work by IMarcus et al.l 
(|2010l ) defines the theoretically anticipated lower enve- 
lope of that distribution. This Letter attempts to sim- 
ilarly define an upper envelope, if one does occur, that 
might separate water planets with no hydrogen/helium 
envelopes from mini-Neptunes with hydrogen/helium en- 
velopes around an icy/rocky or just rocky core. We follow 
the same premise, that well-defined initial bulk composi- 
tions can be modified by late-time giant impacts to reach 
either of the envelopes of the mass-radius distribution of 
super-Earths. 

2. METHOD 

We performed a series of more than 100 simulations 
using the smoothed particle hydrodynamics (SPH) code 
GADGET (jSpringell 12 0051) . modified to rea d tabulated 
equations of state (EOS. IMarcus et~a l. 2009). For more 
details on usin g this scheme f or giant pla netary impact 
simulations see IMarcus et al.l (|2009Ll2010( ). 

We considered targets ranging in mass from 0.5 M 9 to 
5 Mq. The mass of the planets was 50% H2O and 50% 
serpentine (Mg3Si 2 5 (0H)4, ~ 11-15% H 2 by mass) 
with a bulk density of 2.2 — 3.6 g cm -3 . The higher 
bulk densities for the largest mass planets reflect the for- 
mation of dense high-pressure polymorphs of ice and the 
corresponding smaller total planetary radii. This compo- 
siti on is co nsistent with the protosolar ice/rock ratio of ~ 
2.5 (Guillot 1999) and more ice-rich than the composition 
of the dwarf planets and large KBOs in the outer solar 
system ()McKinnon et al.ll2008D . T he serpentine EOS was 
tabulated from the code ANEOS ([Thompson fe Lauson 

mwith input parameters taken from iBrookshaw 
. The ice E OS was the tabulated 5- Phase equation 
of state for H 2 (Scnf t fe Stewart! 12008). Bodies were 



initial ized with fixed temperature profiles from iFu et al.1 
(2010) and all materials are hydrodynamic. Each body 
had 25000-100000 SPH particles. A subset of the sim- 
ulations was run at twice this resolution, with identical 
results. 

The projectile-to-target mass ratios were 1:4, 1:2, 
and 1:1. Because the outcome of giant impacts has 
been shown to depend very strongly on imp act angle 
(jAgnor fc Asphaud [200l IMarcus et all 120091 ). we con- 
sidered impact parameters of 6=0 (0°, head-on), 6=0.5 
(30°), and 6=0.707 (45°). 



3. RESULTS 

3.1. Catastrophic Disruption Criteria and Icy Mantle 
Stripping 

Giant collisions that are near or directly head-on (zero 
impact parameter 6) of sufficient energy lead to disrup- 
tion and net erosion. The criterion for catastrophic dis- 
ruption is an impact condition that leads to a largest 
remnant with half the original combined colliding mass. 
In the disruption regime, the mass of the largest impact 
remnant is a function of the reduced mass kinetic energy 
divided by the total colliding mass, scaled by the catas- 
trophic disruption criteria Qr D , a scaling factor that de- 
pends primarily on the total colliding mass and the im- 
pact velocity. The reduced mass kinetic energy, scaled to 
the total colliding mass is Qr = ^fiV^/M to t, where \i is 
the reduced mass, V\ is the impact velocity, and M to t is 
the total mass. The size- and velocity-dependent catas- 
trophic disruption criteria for weak bodies in the gravity- 
domi nated regime are given by (jStewart fc Leinhardtl 
120091 Equation 2) 



) RD = 10-^fvf 8 . 



(1) 

Here Rci is the radius of a spherical body with mass Mtot 
and a mean density of 1 g cm~ 3 . These disruption crite- 
ria, originally fit to TV-body i mpact simulat i ons b etween 
planetesimals, were shown by lMarcus et al.l (J20Q9) to ex- 
tend into the super-Earth mass regime for bodies of ter- 
restrial composition. In Figure [TJ the fitted Qr D values 
for collisions between icy planets ar e shown as a function 
of Rg i . The disruption criteria of IStewart fc Leinhardtl 
(2009) fit the data for ice-rock planets as well, predicting 
Qrd to within a factor of 2 of the hydrocode calcula- 
tions for head-on (6=0) impacts. The predicted Q RD 
is smaller than the hydrocode calculations for head-on 
impacts; however, the hydrocode calculations do not in- 
clude material stre ngth which biases the valu e of Q RD 
downward (e.g., see ILeinhardt fc Stewart 2009). Oblique 
impacts are discussed below. Thus, these disruption cri- 
teria are largely independent of the internal compositions 
o f the colliding bodies. 

IStewart fc Leinhardtl (|2009l ) found a linear relationship 
between Qr and the mass of the largest remnant scaled 
to the total colliding mass: 

M lT /M tot = -0.5(Qh/Qkd - 1) + 0.5 . (2) 

In FigureEl the mass of the largest remnant, scaled to the 
total colliding mass, is shown as a function of the reduced 
mass kinetic energy. The linear relationship still holds 
for super-Earth planets consisting largely of ice (dotted 
line). Thus, this scaling relationship should be valid for 
all solid planetary bodies, regardless of the details of the 
internal structure. 

In these disruptive simulations, some of the icy mantle 
is removed, leading to an increase in bulk density. The 
rocky core mass fraction of the largest impact remnant 
is given by a power law of Qr/Qr D , 



M COIC /M h . = 0.5 + 0.25(Q R /Q^ D ) 1 ' 2 , 
shown by the dashed line in Figure [2] 



(3) 



3.2. Effect of the impact parameter: Disruption versus 
Hit and Run 



High velocity head-on impacts fall in the disruption 
regime discussed above. As the impact parameter in- 
creases and the impacts become more grazing, a differ- 
ent regime emerges. The hit-and-run imp act events de- 
scribe d by lAgnor fc Asphaual ([2004) and IMarcus et al.l 
(2009) also occur for the planets presently considered. 
This regime is defined by a sharp discontinuity between 
merging impacts at low velocity and high velocity colli- 
sions in which the impactor grazes the target and both 
bodies escape the encounter largely intact. 

Figure [3] shows the accretion efficiency, £, as a function 
of the impact velocity scaled to the mutual escape veloc- 
ity ( Wc 2 sc = 2G(M targ + Af proj )/(i? targ + flproj)). The ac- 
cretion efficiency is defined as the difference between the 
mass of the largest impact remnant a nd the target m ass, 
divided by the mass of the impactor ( Asphaug 2009) : 



_ M\ Y — Aftarg 



M, 



(4) 



proj 



The impact velocity of the transition between accretion 
and hit-and-run depends on the impact angle and the 
projectile-to-target mass ratio, but always falls in the 
range Vi ~ 1.1 — 1.5w csc . For the icy planets, all of the 
6=0.7 simulations transitioned from merging to hit-and- 
run and remain in the hit and run regime for the range of 
impact velocities considered here. The simulations with 
6=0 smoothly transitioned from merging to disruption. 
For 6=0.5, there is a small hit-and-run regime between 
merging and disruption. 

In the hit-and-run regime, catastrophic disruption may 
still occur at extremely high impact velocities. Figure [T] 
includes the calculated catastrophic disruption criteria 
versus combined spherical radius for 30° impacts with 
mass ratios of 1:2 and 1:4; the predicted Q^ D using Equa- 
tion (JTJ are only 20% larger on average when account- 
ing for the different critical impact velocities. The 30° 
impacts with mass ratio of 1:1 and all the 45° impacts 
resulted in hit and run up to 414 sc , and we could not 
derive a <9r D - 

30° impact conditions typically result in <5^ D values 
larger than that for head-on impacts by a factor of ~2. 
However, there is also a non-trivial dependence on the 
mass-ratio of the colliding bodies. The combined factors 
of mass ratio and impact angle are complicated enough 
to warrant a future dedicated study for the boundary 
between the disruption and hit-and-run regimes. 

3.3. Two Models for Mantle Stripping in a 
Differentiated Body 

Combining the results of the simulations of collisions 
between rock/ice super-Earths discussed above with 
the collision s imulations of iron/rock super-Earths from 
IMarcus et al.l (J2OO90 . we construct two general models 
for collisional mantle stripping in a two material, differ- 
entiated planet. In the first one, we begin by calculating 
Qrd from Equation (JTJ, followed by the expected mass 
of the largest remnant from Equation (J2J . In this model, 
we assume that all escaping mass is the lightest material; 
that is, the cores of the two bodies merge. If the com- 
bined mass of the two cores is less than the mass of the 
largest remnant, mantle material is added until the total 
mass reaches M\ T . Given the core masses of the target 
and projectile, M core]ta rg and Af corC!pro j, the mass of the 



final post-impact core is then 
M core = min(Mi r , M cor 



M. 



core, pro 



j)- 



(5) 



In the second model, the first step is again to calcu- 
late Q^ D and the expected mass of the largest remnant 
from Equations (JTJ) and (J2J, respectively. In the accre- 
tion regime, where the mass of the largest remnant is 
larger than the target mass, we assume that core mate- 
rial is accreted first from the impactor followed by mantle 
material. The mass of the final post-impact core is then 



i.targ + min(M c 



j,M lr -M targ ). (6) 



In the disruption regime, where the mass of the largest 
remnant is smaller than the target mass, we assume that 
the material is stripped away from the target in spherical 
shells, beginning with mantle material only. The core 
mass of the largest remnant is 



M„ 



min(M, 



corc,targ? 



M lr ) 



(7) 



A comparison in core mass fraction between these two 
simple models of accretion/stripping and the simulation 
results is shown in Figure |H for both ice- rock bodies 
(a) and rock-iron bodies (b) (|Marcus et al J I2009Q . The 
dashed curve is the prediction from the first model, 
the dash-dotted curve is the prediction from the sec- 
ond model, and the points are simulation results. The 
first model generally overpredicts the core mass fraction, 
while the second model generally results in an under- 
prediction. Given the assumptions of the two models, 
this result is unsurprising. The first model assumes that 
both cores are incorporated into the final body, even in 
the disruption regime. The second model ignores any 
loss of material from the target mantle in the accretion 
regime and that assumes material is lost in shells in the 
disruption regime. In reality, some material will be lost 
from the target mantle even in the accretion regime, and 
in the disruption regime, material is preferentially lost 
from the side of the target suffering the impact, includ- 
ing the core when all mantle material is stripped. 

Despite these discrepancies, the models match the sim- 
ulations fairly well, to within ~15% for ice-rock bodies 
and ~10% for rock-iron bodies. Note that the calcula- 
tions for ice-rock planets were checked for sensitivity to 
the number of particles in the simulation. 

These simple models allow for conservative estimates 
of the limiting case of maximizing the total water con- 
tent of the largest remnant after a giant impact. While 
Figure 0] shows the results for head-on impacts only, the 
model holds also for off-axis impacts in the merging and 
disruption regimes. Hit-and-run impacts result in very 
little material exchange and only minor collisional ero- 
sion of the target and projectile mantles. 

4. DISCUSSION AND CONCLUSIONS 

In previous work, we de monstrated that the catas- 
trophic disruption criteria of lStewart k, Leinna rdt (2009) 
extend into the super-Earth mass regime for rock-iron 
differentiated bodies. Here, we show that the disrup- 
tion criteria still hold for bodies of ice-rock composition. 
Thus, the catastrophic disruption criteria are largely in- 
dependent of internal composition of the colliding bodies, 
holding for gravitationally bound rubble piles as well as 
for icy super-Earths. The catastrophic disruption cri- 
teria (Equation (JTJ) and universal law for the mass of 



the largest remnant (Equation ([5])) are a powerful tool 
for determining collisional outcomes that can easily be 
incorporated into ./V-body planet formation models. 

In the disruption regime, we derive a scaling law for 
the mass fraction of the rocky core following collisions be- 
tween these icy planets. While this relation, as well as the 
similar law for iron mass fraction f ollowing collisions be - 
tween rock-iron planets derived in lMarcus et al.l ((2009), 
is interesting, we have extended the usefulness of such re- 
lationships by showing that a more general law holds for 
how the composition changes in giant impacts between 
differentiated bodies. In our simple model, collisions that 
lead to accretion initially incorporate the heavier mate- 
rials from the impactor while erosive collisions remove 
the lightest material from the target first. This simple 
picture agrees to within ~40% with the simulations for 
rock- iron bodies and is accurate to ^15% for ice-rock 
bodies. This result will be very useful to TV-body planet 
formation model s that track the composition of grow- 
ing p lanets (e.g.. iRavmond et al.ll2006t iThommes et al.l 
[20081) 

Further, our cur rent results, as well as those in 
iMarcus et al.1 (|2009f ) . show that giant impacts between 
bodies of similar composition will never serve to decrease 
the bulk density of the target planet. If the commonly as- 
sumed maximum ice fraction o f 75 wt% for bodies for m- 
ing beyond the snow line (e.g., Mordasini et al. 2009) is 
correct, giant impacts cannot serve as a mechanism for 
increasing the ice fraction. Such impacts would cither 
accrete heavier materials in the same proportion as in 
the target, leaving the ice fraction unchanged, or would 
strip away material from the target's icy mantle, decreas- 
ing the ice fraction. This result confirms scenario (3) for 
the initial states of super-Earth planets bv lValencia et all 
(2007, impacts can cither strip the mantle or provide late 
water delivery in the proportions of the projectile) and 
a maximum radius for water super-Earths in the mass- 
radius diagram. 

One scenario that would seemingly violate this rule 
and result in very underdense planets wo uld be a coll i- 
sion similar to the Moon- forming impact (jCanupl [20041 ). 
followed by a dynamical event that led to the satellite 
becoming a planet, unbound from its parent body. This 
planet would then be highly depleted in core material. 
However, such satellites forming from impact-created 
disks are typically not m ore than a few per cent of the 
mass of the parent body (jCanup et alJl200iT ) and would 
thus not fall in the super-Earth mass regime. Massive 
secondary bodies survive only in hit-and-run events, in 
which case there is only minor collisional erosion of the 
target and projectile mantles. 

Therefore, when considering the possible internal com- 
positions for transiting sup er-Earths with relatively la rge 
radii, such as GJ1214b (jCharbonneau et aD l2009f ). a 
composition of about 75 wt% ice is the highest ice frac- 
tion that is consistent with our current understanding 
of planet formation. If the planet's radius is larger than 
this, the presence of a gaseous hydrogen-helium envelope 
should be assumed. In a sense, the curve on the mass- 
radius relationship corresponding to 75 wt% ice is then 
a likely upper-bound on the radii for terrestrial planets. 
Above this curve, planets should be similar to the ice 
giants of the solar system. 

Not surprisingly, the impact conditions necessary to re- 



move a substantial fraction of an icy mantle are much less 
extreme than those necessary to remove silicate mantles 
IjMarcus et al.l [2010h . While removing an entire rocky 
mantle in a giant impact is unlikely because of the large 
impact velocities necessary, the impact velocities needed 
to strip much of an icy mantle are attainable during 
planet formatio n. Applying the results of this Letter and 
previous work ((Marcus et al.ll2009l [20101. future planet 
formation models incorporating the scaling laws for gi- 
ant impacts can accurately track compositional changes 
during planet formati on. Work analogous to that of 
iMordasini et al.l (J200 9 1 ) can then provide us with statis- 
tical likelihoods for such icy mantle-removing impacts 
during late stage giant impacts. 

The NASA Kepler mission should uncover an unprece- 
dented large number of transiting su per-Earth planets 
((Borucki fc for the Kepler Team! l20lo( ) . We expect that 
many of them will have masses M p < 5Me and rela- 
tively short orbits, ensuring a negligible fraction of plan- 
ets with H/He envelopes. If this expectation turns out 
to be true, the Kepler data should show a clear upper 
envelope in the mass-radius distribution. Otherwise, at- 
mospheric spectro scopy will be needed to distinguish the 
"mini-Neptunes" (Mi ller-Ricci et al.ll2009h . 

The simulations in this Letter were run on the Odyssey 
cluster supported by the Harvard FAS Research Comput- 
ing Group. 
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Fig. 1. — Catastrophic disruption criteria Qp D vs. i?ci (Equa- 
tion {TJ for V\ =10, 20, 40, 60 km s _1 ) and hydrocode results (trian- 
gles: impact parameter b = 0; squares: b = 0.5). Inset shows both 
strength- and gravity-dominated regimes for the size and velocity- 
dependent Qr D - 
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Fig. 2. — Largest remnant mass and core mass fraction vs. scaled 
impact energy for head on impacts. The dotted line is Equation 
J2J. The dashed line is Equation {3}. 
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